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Abstract 

We present a microscopic derivation of the resistive transition equation for 
mixed state of superconductors. This form fits the experimental data of MgB2 
with parameters in agreement with the prediction of BCS superconductivity. 
It also fits the experimental data of high quality untwinned YBCO single 
crystal but with parameters somewhat different from the BCS prediction. A 
discussion in connection with the problem of cuprate superconductivity is 
given. 
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I. INTRODUCTION 

The discovery of high-temperature superconductors (HTS) revived the interest in the 
mixed state physics. One of the most interesting issues is the resistive transition behavior of a 
superconductor. The heuristic pioneering work of Tihkham |IJ showed that by using a proper 
activation energy for flux flow the experimentally measured early field-depending broadening 
of HTS can be accounted for nicely by the model of thermally activated flux flow model 
(TAFF) |§. However, subsequent measurements with higher sensitivity showed that the fit 
was quantitative only down to one tenth of the normal resistance R ~ Q.lR n . A number 
of authors suggested that the much more rapid drop of resistance than the exponential 
predicted by the TAFF model might be due to a freezing / melting transition between a 
vortex liquid and a vortex solid of some sort ||]. Koch et al. H and subsequently Gammel 
et al. || carefully tested experimentally the predictions of the 'vortex-glass model' proposed 
by M. P. A. Fisher by measuring the I ~ V curves and resistance of YBCO epitaxial 
thin films and single crystals with high sensitivity. They found that for each filed B the 
isothermal I ~ V curve shows a power-law behavior at a uniquely defined temperature 
T g (B) as 

V oc I n . (la) 

At T > T g (B) , the curvature on a In V vs In / plot is positive corresponding to the form of 
TAFF models prediction 

V oc sinh fj^j . (lb) 

In contrast, at T < T g (B) one finds a negative curvature with the characteristic predicted 
by vortex-glass model 0] 

(lc) 

Since T g (B) is defined as the dividing point between temperatures for which the linear 
resistance (i.e., R in the limit of / — > 0) is zero, as indicated by Eq. (lc), and that in 

2 



V oc exp 



Jt 
J 



which it is not zero, as Eq. (lb), this transition temperature T g (B) is operationally very 
much the same as the melting temperature of vortex solid T m (B) measured in (nearly) ideal 
crystals. Kwok et al. observed a sharp 'kink' or 'knee' in the magnetic-field-broadened 
resistive transition at R/R n < 0.12 in a high-quality untwinned YBCO single crystal. This 
behavior obeys the angular dependence expected from the Lindemann criterion of vortex 
lattice melting ||. The resistive transition broadening as well as the 'irreversibility line' 
were also observed in the recently discovered MgB2 with T c ~ AOK |§, though it shows 
typical behavior of a phonon-mediated BCS superconductor through the B isotope effect 
0- 

Recently, it is found that the resistive transition of both HTS and MgB 2 can be well 
described with a equation of the normalized form 



R 



— = exp 

tin 



(2) 



^-nii + Vi-Xi^eiTi-T) 

i=l 

with Xi, yi, 7i and Tj the normalized current, voltage, symmetry-breaking factor and critical 
temperature respectively defined as 



Id(T, B) 



x 2 = 



I c0 (T, B) 



R 

yi = 



7! = In 



Id(T, B) 
R n (T, B) 



R 



y ^R } 



Rf(T, B) 
T x = T c (B) , 
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Ueo(T, B) 
U c (T, B) 



kT 

T 2 = T m (B) 



(3) 



where Id is the depairing current and I c q is the critical current of vortex solid for overcoming 
the activation energy barrier U c (T, B). 9 (x) is the Heaviside function and Pi are exponents 
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In present work we try to show the connection of the resistive transition equation (2) 
with the basic Ginzburg-Landau (G-L) theory |12|which can be derived from the microscopic 



theory of inhomogeneous superconductor as shown by Gorkov ||13|| . 

In sections II and III we shall derive the two terms in the right hand side braket of Eq. (2), 
i.e., i = 1 and 2, in connection with the normal-superconducting state (N-S) transition and 



flux pinning respectively. The derived equation will be compared with the experimental data 
of MgB 2 and YBCO single crystal in section IV. A discussion about possible implications 
is given in section V. 

II. TRANSITION NEAR B c2 (T) 

The basic frame for describing the physics of superconducting mixed state is the 
Ginzburg-Landau (GL) free energy density / in the form Jl2| 



/ = / n o + «|vl>| 2 + ^|^| 4 + ' 



% c 



2 h 2 



2m* 

where f n o is the free energy density of normal state, \& (r) is the complex order parameter, 
A is the vector potential. Gorkov found that the GL theory based on Eq. (4) is derivable as 
a rigorous limiting case of the BCS microscopic theory with $ (r) proportional to the local 



value of the gap parameter A (r) and the effective charge e* in Eq. (4) equal to 2e [13 



Within the framework of GL theory Bardeen and Stephen studied the flux flow in a pinning 
free mixed state and obtained the flow resistance pf as [[14J 

p f _2^B faV B ^ B 

Pn $ \U B c2 ~ B c2 > y } 

where a is the radius of normal core in the local model, nearly equal to the G-L coherence 
length £ which can be expressed in terms of Gorkov's derivation with BCS superconductivity 
as [|13[ 



? 1 ; 2m* \a(T)\ I2^m* \ TJ ' K ' 

with ( (3) the Rieman zeta function and (3 C = 1//cbT c . 

For phonon mediated BCS superconductors the energy gap A (T) near T c is shown as 

^'-iifi-I), P) 



A J V T c 

with Ao = A (T = 0) = 1.76ksT c [I2,IS|. Further study on the dependence of A on field 
and Cooper pair velocity v s finds (T5| 



57T 



2 / T 



[arrof ~|^i-^j--(PM 2 ^ (8) 

[&A (T,0)] 2 -I^Mf) 2 v 2 s . 



Inserting (6), (7), (8) into (5) and considering the identity 

In (1 — x) — — — > n—1,2,..., with - 1 < x ^ 1, 

We find 

A 2 (T,B,0) 



(9) 



p f (T,B,v s ) =p n (T,B)e 71 exp In 
« p n (T, 5) exp <^ -71 



A 2 (T,B,v 

2 (hk F f v 



2 2 

5 



3 7o A 2 (T,B,0) 

Using the relation J s = en s v s with the density of superconducting electrons n 
2n (1 - T/T c ) « 2nA 2 (T) /ttA 2 ,, one gets Eq. (10) in the form 



(10) 



Pj (T,B,J) = p n {T,B) exp 



-7,(1-1 



(ID 



with 



•J A 



2 7o e 2 n 2 A 6 (T, 5, 0) 



oc 7 iT c 2 1 



(12) 



(^ F ) 2 7rA^ " \' T c 

Eq. (11) is the same form of Eq. (2) in the pinning-free (72 = 0) case with the evaluation 
of average super current J s as 

E(J) 



J, = J 



(13) 



Pn (T, B) 

III. RESISTANCE DUE TO FLUX CREEP 

Working in the London limit, Nelson showed that according to the GL free energy Eq. 
(4) a system of iV flux lines with a field H along the z direction in a sample length L can 
be described with the free energy represented by the trajectories Irj (z)| of these flux lines 
16Hl8|. Considering further the pinning potential Vp ( r ) arising from inhomogeneities 



and defects in sample [f7[,19j20|, the free energy of such a sample with iV flux lines can be 
expressed as 



1 N 

1 ;-!•"> 



d Tj (z) 



dz 



1 r L N 

2 m Jo 3=1 







[ L v P 
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dz. 



(14) 
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Here V (r^) = V y r j — r 3 - J = 2e K (rij/X ab ) is the interaction potential between lines 
with in-plane London penetration depth X a t and Kq (x) is the modified Bessel function 
Kq(x) ~ (n/2x) 1 ^ 2 e~ x . E\ is the linear tension of flux line and e ~ ($ /4vrA a fc) 2 is the 
energy scale for the interaction. 

Thermally activated flux motion is considered as the sequence of thermally activated 
jumps of the vortex segments or vortex bundles between the metastable states generated by 
disorder. Every elementary jump is viewed as the nuclearation of the vortex loop, and the 
mean velocity of the vortex system is determined by the nuclearation rate R,[19|j2T 



v oc exp (SF/kT) 



(15) 



Here SF is the free energy for the formation of the critical size loop or nucleus which can be 
found by means of the standard variational procedure from the free energy functional due 
to the in-plane displacement u (z) of the moving vortex during loop formation 

2 



loop 



U 



dz 



1 

-£/ 
2 



d u (z) 



dz 



+ V P 



u (z) - (f L + fr,)-u 



(16) 



where fi=Jxe z /cis the Lorentz force due to applied current J and f v is the viscous 
drag force on vortex, f v = —r]V vortex , with v vortex = d u jdt and viscous drag coefficient 
r\ ~ ($o-£>c2) / (PnC 2 ) as estimated by Bardeen and Stephen fL4| . 

Equation (16) is similar to the basic equations of collective pinning model and vortex- 



glass as well as Bose-glass models |21] for considering the vortex dynamics at low tempera- 
tures. The only difference is the omission of f v in the latter, which is quite resonable since 
the velocity of vortex v vorte x usually is very small in the case of low temperature flux creep, 
where one finds f v fi. However, we have to consider the term f v in our equation (16) not 
only because the fact that viscous drag force does ever accompany the vortex displacement 
as shown experimentally by Kunchur et al. p2 but also for the need in our task to describe 
the current-voltage characteristic of vortex system more accurately for the case across the 
irreversibility line where fi and f v are comparable in magnitude, so now we have 
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loop 



U 



dz 



2 



d u (z) 



dz 



+ V P 



U (z) - f s - U 



(16') 



with 



fa = h + Jr,= x e 2 



and 



Pf 



(17) 



where we used the relation E =v x B and Eq. (5). 



The free energy functional F loop 



u 



contains two parts. The last term in the right 
hand side of Eq. (16) represents the net energy gained from the applied current which is 
proportional to J p in Eq. (17) and the rest are corresponding to the increase of the elastic 
energy during the loop formation. As argued by Fisher et al. 0, in a bulk superconductor 
vortex lines are extended one- dimensional (ID) objects and the response of vortices to an 
applied current can be described by vortex loop formation of area S = L± ■ L z ~ L\ and 
elastic energy increase ~ U c ■ L 9 ± , where L± is the transverse displacement of vortex-line 
segment of length L z , with scaling relation L±_ ~ L\. 

Consequently, the free energy of loop formation can be estimated in the dependence of 
L i as 



Fioop (L±) ~ U c Lj_ — Jp—L 1 ^. 



(18) 



For a given value of J p , F\ oov (L±) first increases and then decreases with increasing L± (since 
k > 9 and J p < J c ) . The barrier energy is the maximum value of loop formation free energy 



5F=U (J p ) = Fi oop [L* ± (J p )] 



with definition of L* ± (J p ) 
dFi oop (Lj 



dL 



0, L* ± (J p 



L i =L* 



Combining Eqs. (18), (19) and (20) one finds 



cBU c 



V(*-«) 



(19) 



(20) 



(21) 



where 



eu c c 



and n 



e 



2C-1 



(22) 



/t — 9 2-C 

The current dependence (21) of the barrier U to current-driven flux motion implies a current 
voltage characteristic of the form 



E (J) = pjJexp 



U c ( J c 



kT V J, 



(23) 



This form in one hand turns to the pinning-free flux flow resistive regime as J, J p ^> J c and 
clearly goes to zero in the other hand as J p ~ J — > 0, with no linear resistance term. 

In real samples, the critical size of loop formation L* ± (J p ), and thus [/ (J p ) may always 
manifest a nonzero resistance at sufficiently small measuring current with high enough accu- 
racy of measurements. Considering this real size effect one finds a general normalized form 
of the current-voltage characteristic in the form (see Appendix) 



y = xexp [-7(1 +y-x) p ], 



(24) 



with 



U r (J, 



7 kT \J L 



J E(J) 

x = -r, y = — 
Jl PfJL 



p = p. 



(25) 



Where Jl is the transport current density corresponding to the case where the critical size 
of loop formation is equal to the sample size L as 



c9U c 



Ji 



E{J L 
Pf 



Jr 



Jl 



E{J L 
PS 



^ = L K - 
Jl 



(26) 



Together with the former equation (11) in section II, equation (24) is just the normalized 
form of the resistive transition equation (2) in the case of temperatures below the melting 
point of vortex solid T < T 2 = T m (B) with p 2 = p and J c0 (T, B) = (T, B). 

IV. COMPARISON WITH EXPERIMENTS 

In this section we compare the above derived resistive transition equation of mixed state 
in the form of Eq. (2) with the experimental data of several typical superconductors with 
high T c . 



A. MgB 2 

The recently discovered superconducting MgB 2 has a high critical temperature T c ~ 4QK 
comparable to the cuprate LaSrCuO . Meanwhile, it shows isotope effect of phonon- 
mediated superconductivity. Finnemore et al. measured the transport and magnetic prop- 
erties of sintered pellet of MgB 2 and find the Ginzburg-Landau parameter k ~ 26 0. In 
Fig. 1 we show the comparison of the resistive transition equation (2) with their experimen- 
tal data of the temperature dependent resistance of MgB 2 from 300K to 1.9K in different 
applied fields. 

B. YBa 2 Cu 3 7 - 5 (YBCO) 

It is widely believed that the mechanism of superconductivity in high-T c cuprates may 
be essentially different, from the familiar s-wave type pairing on which conventional BCS 



theory is based ||22|| . Kwok et al. studied the width and shape of the resistive transition 
of untwinned and twinned single crystals of YBCO in fields up to 8T || .This is one of 
the first evidences which came from transport measurements for vortex melting in YBCO. 
The samples are of high quality and near optimum doping confirmed by their zero-field 
resistive transition of R zero > 92.0K and AT C (10% — 90%) < 0.2K. A "knee" in R(T) 
curve is clearly seen. We compare our resistive equation (2) with their experimental data of 
untwinned YBCO crystal in Fig. 2. 

Agreement is fair in both cases. 

V. DISCUSSION 

Though the resistive transition equation (2) fits the experimental data of both MgB 2 and 
YBCO crystal, there are still essential differences in the parameters entered into Eq. (2) for 
these two kinds of materials. In Fig. 1 we see the parameters used to fit the experimental 
data of MgB 2 by Finnemore et al. || follow the relations Id oc 7°' 5 [T c (B) — T] ' and 
7a = In [R n (T, B) /R f (T, B,J ^ 0)] = In [B c2 (T) / B] in agreement with Eq. (17) and the 
flux-flow resistance equation derived by Bardeen and Stephen Q for BCS superconductivity. 
In contrast, the parameters to fit the experimental data of YBCO crystal by Kwok et al. 
§show the relations I d oc 7?- 5 [1 - T/T c (B)f 3 and 71 = In [R n (T, B) / R f (T, B, J -> 0)] = 
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In [B C 2 (T) /B} m , with m = 0.845 -0 ' 59 — 0.07 which are somewhat different from Eq. (17) and 
the results of Bardeen and Stephen in the exponents. At present a correct microscopic 
theory for the high-T c cuprates is still a challenging problem in condensed matter physics. 
While it seems that the established superconductivity of cuperates is of d x 2_ y 2 symmetry, 
the next question is whether it can be described by a BCS- like theory suitably modified 
to in elude a d-wave gap or it is of some non-BCS origins ||23|| . On the one hand, there 



are recent reports on the evidences for strong electron-phonon coupling in high-T c cuprates 
by the analysis of photoelectron spectra [24| and the unconventional isotope effect |25| 



as well as the identification of the bulk pairing symmetry in high-T c hole-doped cuprates 
as extended s-wave symmetry with eight-line nodes and as anisotropic s-wave in electron- 
doped cuprates [26| which seem favouring the BCS-like machanism. On the other hand, 



the well established unusual properties of both the normal and superconducting states of 
cuprates seem related to the fact that the cuprates are doped Mott insulators and stimulate 
some ideas based on new exitations such as spin-charge separation p7|j28[1 , stripes [29|, new 



symmetry relating superconductivity and magnetism [RO] or quantum critical points pTT 



which suggest a non-BCS state. The interplay between theory and experiment promises 
to be mutually benificial, in the best traditions of physics research. In this paper we have 
compared the derived equation (2) only with the data of optimally doped YBCO single 
crystal ||. However, the unusual properties of cuprates appear even more striking in the 
underdoped region. A thorough study on the resistive transition of underdoped cuprates 
may provide new important insights into the nature of cuprates. 
VI. SUMMARY 

Starting from the Ginzburg-Landau functional we derived a resistive transition equation 
for the mixed state of superconductors. This equation in its general form agrees with the 
experimental data of superconducting MgB 2 pellet and optially doped untwinned YBCO 
single crystal but the parameters to fit the data of cuprate are somewhat different from the 
prediction based on BCS theory. 
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APPENDIX 

In this appendix we provide the derivation of Eq. (24). Substituting Eq. (17) into the 
bracket on the right-hand side of Eq. (23) and taking its logarithm we get 



Jr. 



-i/n 



(A.l) 



with Jf = E (J) /pf. Since the critical size of loop formation L* ± is limited by the sample 
size L as L* ± L, one finds always 



U{J)^U{J L ) 



with the definition 



c9U c 



Jl~ 



E{J L 
Pf 



(A.2) 



(A.3) 



Here we have used the relations in Eq. (17) and Eq. (20). Thus, Eq. (A.l) can be expressed 
in the form 



J- J 



f 



U, 



x/n 



hi 



J J. 



J, 



1 -i/m 



(A.4) 



where Jif = E (Jl) /pf, which is much smaller than Jl and 

- In (J^) 1 ^- In (J/J^l 
h = —± 



|*|<1. 



(A.5) 



In {Jl/Jlj) 1 '" 

Using the approximation (1 + h)^ 1 for \h\ < 1, finally we find Eq. (A.l) in the form 



x — y = 1 — In — 

\y 



X 



1/p 



7 



-i/p 



(A.6) 



which is exactly the current- volt age characteristic Eq. (24) with the definitions 
Jl fU r \ ( X ^ /tt \ / t\p 



7 = In 



x = 



J, 

Uc_ 

kT 



Lf 



kT J \ Jr - J, 



-1/p 



-1/p 



111 



Lf 

J L \ 1/p fJ 



Jlf ) V Jr 

1/p 



Uc_\ Jc 

kT) \J L/ 

1 J J 

2 Jl — Jlj Jl 



E{J) 



'' 1 kT) \ n J L fJ V JoPf J 2 (J L - J Lf ) p f ~ J LPf 



(A.7) 
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FIGURES 

FIG. 1. Comparison of Eq.(19) with the experimental resistance data of MgE>2 samples mea- 
sured by Finnemore et al. ||. (a) A full view in the applied field of 9T. (b)Resistive transi- 
tions in different applied fields. □ O Av* O +"A"X denote experimental data and lines denote 
the theoretical curves of Eq.(19) with corresponding applied fields. The parameters in Eq.(19): 
1 U(T,B) oc [T m (B)-T]°- 8 (B + 5A6)- A -\ 2 I c0 oc T m (B) — T; I d oc [T c (B) - T} 1 ' 5 ^ 5 . 



3 T m (B) = T c (0) 1- (5/21.7) 



x0.84 



TAB) 



T c (0) (1 - 5/22.4), T c (0) = 40.2K. 4 



70 = In [Rn (T, B) /R f (T, B,J^ 0)] = In [B c2 (T) /B), where B c2 (T) = 0.6 [T c (0) - T] 



0.98 



FIG. 2. Comparison of Eq.(19) with the experimental resistive transition data of un- 
twinned YBCO crystal measured in different applied fields for 5// c by Kwok et al. 
||. DO^v** +"A _X denote experimental data and lines denote the theoreti- 

cal curves of Eq.(19) with corresponding applied fields. The parameters in Eq.(19): 1 
U(T,B) cx [T m (B)-T] 3 B- 4 - 22 . 2 / c0 oc T m (5) — T; I d oc T c (B) [1 — T/T c (5)] 2 ' 3 7q- 5 . 3 
T m (B) = T c (0) [l - (5/1200) a36 l, TAB) = T c (0) [l - (5/3.6 x 10 6 ) 03 ] , T c (0) = 97K. 4 
70 = In [Rn (T, B) /Rf (T, 5, J — > 0)] = In [5 c2 (T) /5] m , where m = 0.845- 059 -0.07 and 
5 c2 (T) = [T c (0) - T] 3 ' 3 . 
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